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ABSTRACT 



Context. Numerous surveys are currently underway or planned that aim to exploit the strengths of the Lyman-alpha emission line for 
cosmological purposes. Today, narrowband imaging surveys are frequently used as a probe of the distant universe. 

Aims. To investigate the reliability of the results of such high-z Lya studies, and the validity of the conclusions that are based upon them. To 
determine whether reliable Lya fluxes (F Lya ) and equivalent widths (W Ly „) can be estimated from narrowband imaging surveys and whether 
■ any observational biases may be present. 
Methods. We have developed software to simulate the observed line and continuum properties of synthetic Lya galaxies in the distant universe 
by adopting various typical observational survey techniques. This was used to investigate how detected F Lya and Wi# a vary with properties 
of the host galaxy or intergalactic medium: internal dust reddening; intervening Lya absorption systems; the presence of underlying stellar 
populations. 

Results. None of the techniques studied are greatly susceptible to underlying stellar populations or the relative contribution of nebular gas. We 
Q ' find that techniques that use one off-line filter on the red side of Lyo- result in highly inaccurate measurements of Wiy„ under all tests. Adopting 
, two off-line filters to estimate continuum at Lya is an improvement but is still unreliable when dust extinction is considered. Techniques that 
C/5 employ single narrow- and broad-band filters with the same central wavelength are not susceptible to internal dust, but Lya absorption in 



Or 



c3 



the IGM can cause W^ ya to be overestimated by factors of up to 2: at z = 6, the median W Lya is overestimated by ~ 25%. The most robust 
approach is a SED fitting technique that fits E B _ V and burst-age from synthetic models - broadband observations are needed that sample the 
UV continuum slope, 2175A dust feature, and the 4000A discontinuity. We also notice a redshift-dependent incompleteness that results from 
DLA systems close to the target LAEs, amounting to ~ 10% at z = 6. 

Key words. Methods: observational - Cosmology: observations - Galaxies: high-redshift - Galaxies: starburst 



1. Introduction 

Galaxies hosting young, violent starbursts are expected to ac- 
commodate numerous massive, hot stars that are bright in the 
far ultraviolet (UV) regime. UV photons with wavelengths 
shortwards of the Lyman break ionise their local interstellar 
media (ISM), resulting in recombination nebulae bright in hy- 
drogen lines. Two-thirds of the Ly man-continuum photons are 
reprocessed as Lyman-alpha (Lya) under the assumption of 
case B recombination, and it could be expected that young star- 
burst regions should be consistently Lya-bright. 

Early observations of local star-forming galaxies with the 
International Ultraviolet Explorer (IUE) demonstrated this not 
to be the case. In small samples of relatively un-evolved galax- 
ies (low dust and metal content) Lya was frequently shown to 
be absent or unexpectedly weak (eg. Meier & Terlevich 1982 
Hartmann et al. 119881 Calzetti & Kinnev [T992l l. Although dust 



could be invoked to explain the absence of this FUV emission 
line, Giavalisco et al. ( 1996) demonstrated that pure dust ex- 
tinction could not explain the weak Lya fluxes from any of 
the galaxies in the IUE samples. Instead results suggest that 
Lya photons have decoupled from the non-resonant continuum 
radiation, resulting in line-photons escaping the host over sig- 
nificantly extended path-lengths; thereby increasing their sen- 
sitivity to dust (Neufeld Q990J. On the other hand, Neufeld 
( 1991 ) also found that resonance scattering can cause less at- 
tenuation of Lya photons compared to continuum if the ISM 
is multiphase: if dusty, neutral clouds are embedded in a dust- 
free, ionised ISM, resonance scattering can reflect Lya photons 
from the surface of the clouds, allowing them to diffuse out 
of the ISM relatively unattenuated. Radiative transport in such 
multiphase, clumpy media has more recently been studied by 
Hansen & Oh (2005 1 who found Lya equivalent widths {Wi_ya) 
could be 'boosted' by factors of 2 - 3. 
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In the local universe, the observational situation was 
further-complicated by Lya spectroscopic observations using 
the Goddard High-Resolution Spectrograph ( GHRS) and Space 
Telescope Imaging Spectrograph (STIS) both onboard Hubble 
Space Telescope (HST). In a sample of 8 local starbursts, 
Kunth et al. ( 1998 1 found that when Lya is seen in emission, 
blueshifted Lya absorption features are also common, leading 
to P-Cygni-like profiles. In the same study, O i and Si n UV ab- 
sorption lines were also blueshifted with respect to the ionised 
gas, indicating large-scale outflows of the ISM with velocities 
up to 200 km s _1 . Mas-Hesse et al. (|2003l demonstrated that 
varying Lyo- profiles from starbursts are well explained by an 
evolutionary super-shell model; over the lifetime of a starburst, 
absorption, pure emission, and P-Cygni emission phases will 
all be observed, depending on the properties of the outflow and 
viewing angle. If Lya escape is associated with winds and pho- 
tons can diffuse through H i, low surface brightness Lya emis- 
sion may be expected to span large areas, significantly more 
extended than the starburst itself as probed by FUV emission. 
Hence imaging becomes an efficient method by which the line 
can be isolated. This was the motivation for our HST survey 
using the Advanced Camera for Surveys (ACS; Kunth et al. 
2003 1. Such Lya imaging studies are extremely sensitive to 
how the continuum is subtracted and much supporting data is 
required in order to estimate the continuum flux at Lyo-. We 
found in Hayes et al. ( 2005 1 that at least three off-line imaging 
observations are necessary in order to estimate the continuum 
flux at Lya. We determined it necessary to sample the UV con- 
tinuum slope and 4000A break in order to disentangle the ef- 
fects of age and reddening, allowing us to estimate the flux due 
to continuum processes in the on-line filter from synthetic mod- 
els. In our sample, luminous blue compact galaxy ES0338- 
IG04 shows regions of diffuse Lya emission spanning several 
kpc, and patchy structure consistent with an inhomogeneous 
and outflowing ISM. In the diffuse regions W^ya is shown to be 
very large (> 200A) where continuum flux is weak, but inte- 
grated over the whole galaxy W\_ ya is ~ 20A, an order of mag- 
nitude lower. Damped absorption with W\_ ya ~ -40A is also ob- 
served in some central regions, demonstrating the sensitivity of 
Lya to small-scale variations in the ISM. 

Due to its sensitivity to star-forming activity, high natural 
equivalent widths (up to 240A; Chariot & Fall 1993) sensitiv- 
ity to dust, and convenient spectral positioning in the FUV, Lya 
has long been considered a competitive tracer of primeval high- 
Z galaxies as they form their first generation of stars. Partridge 
& Peebles J1967b suggested that in a collapsing Hi halo, as 
much as 7% of the gravitational energy could be radiated in 
the Lya line. While surveys that target continuum emission or 
utilise dropout techniques may have been highly successful in 
recent years, they all share one common weakness: a luminos- 
ity bias in favour of galaxies with strong continuum. Such a 
bias leads to incompleteness at low luminosity and these sur- 
veys would miss the numerous low-mass, star-forming objects; 
the dominant species in the paradigm of hierarchal galaxy for- 
mation. 

The early days of high-z Lya studies were largely unsuc- 
cessful (see Pritchet 1994 for a review) and it was not until 



the late 1990's that searches for high-z Lya-emitters (LAE) be- 
came fruitful. Despite the complications, high-z Lya-emitting 
objects are now routinely being discovered by imaging tech- 
niques (Cowie & Hu fT9%l Rhoads et al. l^UUOl Fynbo et al. 
2003 ; Kodaira et al. l2003l Yamada et al. l2005t . fields surround- 
ing massive objects (eg. Monaco et al. 2005 ), and spectroscopic 
observations (eg. Kurk et al. 2004 1. After much initial disagree- 
ment, the predicted and observed luminosity functions (LF) of 
LAEs at high-z may be converging (Le Delliou et al. [2005 1. 
Interestingly however, a fraction of the LAEs uncovered in the 
high-z universe have very large measured equivalent widths 
(Cowie & Hu ll998l Malhotra & Rhoads 120021 In the case of 
Malhotra & Rhoads (|2002i, their LAE population at z = 4.5 
was shown to have a median Wiy a of 400A, 65% greater than 
the maximum value of 240 A (Chariot & Fall, 1993 1 that can be 
generated by pure star-formation models with normal metallic - 
ities and IMFs. These authors interpreted this to imply that ei- 
ther these objects were type II quasars or star-forming galaxies 
with very top-heavy IMFs or zero-metallicity stars. Subsequent 
X-ray observations with the Chandra satellite ruled out the pos- 
sibility that these objects could be AGN (Wang et al. 2004 1. 
Jimenez & Haiman ( 2006 ) demonstrate how this (and other ef- 
fects) can be explained if 10-30% of the stars in these galax- 
ies are primordial. 'Normal' starburst conditions could still ex- 
plain such high equivalent widths but the system would have to 
be perfectly configured so as to allow Lya photons to leak out 
while dust blocks much of the stellar UV radiation. Another 
possibility is that this is the result of some other astrophysi- 
cal effect, external to the host galaxy that may cause certain 
observational techniques to overestimate W\_ ya when such nar- 
rowband imaging techniques are used. 

While Lya may provide a 'clean' probe of the high-z uni- 
verse, the common observational trade-off is evident: spectro- 
scopic studies may be rich in information regarding the line, 
while narrowband imaging studies may be efficient for de- 
tections. Future generations of high-z-optimised integral field 
units (eg. the Multi Unit Spectroscopic Explorer (MUSE) 
bound for ESO's VLT) have the potential to provide significant 
advances in this area. However, in the cases where narrowband 
imaging alone is used to derive information about the line, it 
is essential to examine exactly how efficient the technique is in 
doing so. To our knowledge, no such comprehensive study has 
previously been performed and this article represents the first 
step in doing so. We here simulate how various survey tech- 
niques may estimate Lya detection properties - line flux (F L yQ .) 
and equivalent width {Wt, ya ) - of galaxies at high-z. 

The paper is organised as follows: in Sect. |2 we describe 
the software; in Sect. [5] we present the simulations we have 
performed; in Sect. |4] we present and discuss some of the more 
important results; and in Sect.[5]we present our concluding re- 
marks. We assume a flat cosmology with Qm = 0.3, Qa = 0.7, 
Ho = 70 km s Mpc~' throughout. 

2. The software 

The Lya galaxy simulation software follows a simple three- 
step prescription: (i) the restframe spectral energy distribution 
(SED) of a the test galaxy is generated and the Lya line is 
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added; (ii) the SED is redshifted and the effects of intervening 
material (IGM absorption) are applied to the spectrum; and (iii) 
the spectrum is convolved with various filter response profiles 
and fluxes are computed. 

2.1. SED generation 

Restframe SEDs of starburst galaxies are generated from the 
Starburst99 (hereafter SB99; Leitherer et al. fl999l Vazquez & 
Leitherer 2004 ) synthetic spectral models. No nebular emission 
lines are thought to be strong enough to significantly contribute 
to fluxes in the UV filters considered here (see Sect. I2.3I for 
a description of the filters) - filters are broad and continuum 
dominated, and in Hayes et al. J2005> we concluded that opti- 
cal emission lines had negligible impact upon our study. Hence 
the SB99 spectra are used unmodified, free of nebular emis- 
sion lines but including nebular continuum emission. The soft- 
ware was designed for flexibility and to cover as wide parame- 
ter space as possible, enabling the user to perform a wide array 
of parameter dependency tests. Spectra can be selected from 
the 1999 SB99 data-release, choosing with burst modes (in- 
stantaneous or continuous), composition (stellar-only or stel- 
lar+nebular), metallicity, initial mass function, and age. In ad- 
dition to the standard set of spectra, SED modification code 
was written to allow the selection of: 

1. Departures from the single stellar population (SSP). It is 
well known from studies of the low- and high-z universe 
that star-forming galaxies often exhibit composite spectra, 
with contributions from the current starburst and an under- 
lying, older stellar population. Our method for testing the 
effects of such population mixing is to allow any single 
SB99 spectrum to act as a template for a secondary pop- 
ulation. This secondary spectrum is renormalised so it con- 
tributes a specified fraction of the spectrum under consid- 
eration at a user-specified wavelength, and is added to the 
spectrum or set of spectra. 

2. Variations in the contribution of nebular gas. By subtract- 
ing out the nebular component from corresponding pairs 
of spectra, the nebular gas spectrum as a function of age 
is recovered. The gas spectrum can then be re-scaled and 
added back to the stellar-only component. Thus a normali- 
sation factor of corresponds to a stellar-only spectrum, 1 
corresponds to the standard stellar+nebular spectrum, and 
factors > 1 imply the gas spectrum has been boosted. 

Restframe spectra are rebinned to 1A steps in wavelength. 
After generation of the template spectrum, the whole SED is 
renormalised to L,\ at 1500 A based upon input star-formation 
rate (S FR) using the relationship 

L^ISOOA) = 9.5238 x lO^SF^Moyr 1 ) (1) 

adapted from Kennicutt ( 119981 1. The all-important Lya line of 
specified equivalent width is added to the spectrum at wave- 
length of 1216A. At the rest wavelength of Lya, 1A corre- 
sponds to the velocity of 250 km s _I . In addition, pure imaging 
studies aren't sensitive to line profiles and our aim is simply 
to investigate how well intrinsic line-strengths are recovered 



(whatever their origin or profile). Hence all the energy is de- 
posited in one 1 A bin. See Sect. l2.3l for a more detailed discus- 
sion of this. 

Lya is the only Lyman series feature added to the SED - the 
O vi-Ly/3-Cn complex between 1032 and 1038A does not fall 
within any of the filter bandpasses considered here (Sect. l2~3l . 
The Civ/11549A absorption feature (although frequently ob- 
served with P-Cygni profiles), common in star-forming galax- 
ies, does fall centrally in one of the bandpasses we use. In a 
sample of 45 local galaxies observed with the IUE, Heckman 
et al. (1998) find a (C iv+Si iv)/2 median equivalent width of 
-4. 6 A, with the absorption features never stronger than -10. 1 A. 
Crowther et al. y006l re-measured these equivalent widths, 
finding an offset of 1-2A relative to the Heckman et al. re- 
sult. The filter bandpass that the C iv feature falls in is broad 
((1 + z) ■ 300A, centred at 1500A; see Sect. I2.3t so the pres- 
ence of the strongest locally observed C iv feature would affect 
the integrated flux by only ~ 3% with no redshift dependence. 
Moreover, these authors found that this equivalent width is pos- 
itively correlated with metallicity. Since strong Lya emission is 
often interpreted as a sign of low -metallicity stars, we decided 
not to apply the C iv feature to our SEDs. 

In order to complete the restframe SED, it can be reddened 
using the Galactic laws of Seaton J1979i . or Cardelli et al. 
GHH3l, the SMC la w fit o f Prevot et al. fT985l . the LMC law 
of Fitzpatrick et al. dl985> . and the Starburst law of Calzetti et 
al. J20001 

2.2. Effects of cosmic distance 

The restframe SED is first shifted to the desired redshift and 
the luminosity distance computed from the formula 

d h = (1 + z)^- f = . = (2) 

where Ho is the Hubble constant and Qm and Q.^ are the en- 
ergy densities of total matter and cosmological constant, re- 
spectively. The luminosity at each wavelength is then converted 
to flux using the inverse square law. 

Given that we want to study all the possible astrophysical 
effects, including possible extreme cases, we chose to treat the 
effects of intervening H 1 clouds by generating random distri- 
butions of individual clouds, not by applying some general av- 
erage prescription (eg. Madau 1995 1. The effect of interven- 
ing H 1 clouds (Lya forest (LAF) to damped Lya (DLA) sys- 
tems) is implemented by first assuming the number density of 
absorbing clouds as a function of redshift takes the form of 
the power-law N(z) = No(l + z) y in the redshift range to z 
where No = 0.07+°;^ and y = 2.45 (Peroux et al. l2003l . 
Firstly the total number of clouds in the given redshift range 
is generated by biasing a uniform variate pseudo-random num- 
ber (PRN) by the distribution function, within the observational 
constraints. For each cloud, a position in redshift space is gen- 
erated in the same manner: by weighting a PRN by the dis- 
tribution function. The column density («hi) of each cloud is 
generated by assuming the distribution obeys the power-law 
N(n m ) <* (Rao et al. l2005l for ft = 1.4 + 0.2 and column 
densities in the range 13 < log(nHi[crrT 2 ]) < 22. 
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Once the redshift-column-density distribution is generated, 
the effect on the SED of each cloud is determined from its 
equivalent width in absorption (W a b s ). We compute restframe 
W a bs from the effective optical depth at line-centre (tq) follow- 
ing the curve of growth method (Spitzer 1978). After cosmo- 
logical scaling (1 + Zcioud), W a b s and the continuum flux-density 
at ( 1 +Zcioud) - /lLyff,0 are then used to remove the requisite amount 
of flux from the appropriate bins in the redshifted spectrum. In 
cases where the absorption is damped {W d bs is greater than the 
bin size), the flux in that bin is set to zero and flux is symmet- 
rically removed from neighbouring bins until the required W a b s 
has been met. 

Significant attention has been paid to the possibility of 
reddening of high-z quasar spectra by DLA systems (eg. Pei, 
Fall & Bechtold 1991). However, recent studies have demon- 
strated the typical mean cumulative reddening from DLA sys- 
tems to be small. For example Ellison, Hall & Lira (2005) 
found mean Eb-v < 0.04 at the 3<x level for SMC-type dust 
in a sample of 14 high-z quasars with DLAs, while from a 
sample of 72 DLA quasars from the Sloan Digital Sky Survey, 
Murphy & Liske (2004) found mean E B _ V < 0.02 at 3cr for 
SMC dust. On the other hand, Wild et al. (2006) find evidence 
for significant reddening in Can and Mg n absorbing systems: 
Eb-v = 0.105 for SMC type dust in the selected strongest ab- 
sorbers {Eb-v = 0.066 for their complete sample). Such ob- 
jects are shown to be similar to DLAs with a number density 
of ~ 20-30% that of DLAs at the same redshift. While we do 
not claim that some high-z LAEs would not be heavily dust 
reddened, we chose for this part of the study not to implement 
dust-reddening by high-Hi column-density systems. 

At z = 0, the redshifted SED can be reddened once more 
using the laws of Cardelli et al. ( 1989) or Seaton ( 1979) to sim- 
ulate Milky-Way reddening. 

2.3. Filters and computation of line fluxes and 
equivalent widths 

For the sake of simplicity (our aim was to study astrophys- 
ical processes, not observational complexities) the filters are 
perfect. Defined only by their central wavelength and effec- 
tive width, they provide 100% transmission inside the pass- 
band and zero otherwise. That is, the filters are as close to 
perfect top-hats as possible, without introducing numerical, 
resolution-dependent inconsistencies at the edges. Narrowband 
filters targeting the Lya line are defined to have a full width of 
2% their central wavelength (A c ) while broadband filters have 
full widths of 20% A c . That is, the width of the broadband fil- 
ters scales with A c in approximately the same way as it does 
in the Johnson-Cousins/Bessell and JHKLM systems. The nar- 
rowband filters scale similarly so as not to introduce effects that 
may result from sampling different spectral regions at different 
redshifts. 

One issue that arises when using narrow filters to isolate an 
emission line is the positioning of the line within the filter pro- 
file. That is, is the line centrally positioned and narrow enough 
to allow maximum transmission of the line flux? Chemistry and 
manufacturing processes dictate that filters (particularly nar- 



rowband) cannot be perfectly rectangular and most frequently 
take the form of bell-curves. Whilst photometric redshifts or 
drop-out techniques may be used to remove interlopers (most 
notably [On]/13726A), photo-z's of individual objects do not 
reach the required levels of accuracy to determine the position 
of Lya within a narrow bandpass (typical photo-z accuracies 
of Az = 0.05(1 + z) cannot determine whether the line falls 
within a 2%A C narrowband filter or not). In addition, photo-z 
methods tend to be reliant upon sampling the Balmer/4000A 
break which is lost from optical multiband datasets at z ~ 1, 
something we discuss in another context in Sect. |4] At high-z 
and without spectroscopic data, the selection of LAE candi- 
dates by detection in the narrowband filter becomes the most 
accurate estimate of redshift. Hence effects that result from the 
line falling in the wing of a filter will always be something that 
has to be considered. Starburst galaxies do not exhibit signif- 
icantly broadened spectral lines: concerning Lya, Matsuda et 
al. (|2006 ) spectroscopically measured the widths of 37 LABs 
at z = 3.1, finding Lya full widths (FWHM) in the range 150- 
1700 km s _1 with a median value of 550 km s _1 . This corre- 
sponds to a maximum Gaussian cr of 3 A at Lya (median value 
of ~ 1A). Thus given the 2%A C narrowband filters used here, 
even the broadest Lya lines will be completely transmitted (2% 
rectangular filters drop to 99% of complete line transmission at 
cr = 4.70A ~ FWHM = 2700 km s" 1 ). The bin size of our 
SEDs is 1 A and, since the median cr of the Matsuda et al. study 
is 1A, we feel safe in depositing all of the Lya flux in one bin 
(see Sect. ED- 

'Observed' fluxes in all filters are computed by convolution 
of the profile and the spectrum, and integration. In order to es- 
timate true line-fluxes (F ^ ya ), it is essential to have a robust es- 
timate of the flux due to continuum processes that is present in 
the on-line filter. This estimate of the continuum flux at 1216A 
is central to Lya studies and here we use various methods to 
scale the nearby continuum flux to Lya (described below). 
With estimates of the line and continuum fluxes, equivalent 
widths, defined as W]_ ya = F^ ya /f cont , are computed, where / cont 
is the continuum flux density at line-centre. To this end, a num- 
ber of different commonly implemented and feasible observa- 
tional techniques are employed from which F^ ya and Wi, ya are 
estimated. 

For on-line observations, a narrowband (2%/l c ) filter is 
used, centred at the observed wavelength of redshifted Lya. 
Broadband (20%A C ) filters are positioned at various wave- 
lengths in the redshifted SED in the UV and optical, sampling 
the restframe SED at 1216, 1500, 2200, 3300, and 4400A. That 
is, filter wavelengths and full widths are defined in the rest- 
frame - the same region of the SED is sampled at all redshifts. 
Using these broadband filters, continuum flux at Lya is esti- 
mated using four techniques numbered (#l-#4) in order of the 
central wavelength of the reddest filter used. A schematic rep- 
resentation can be seen in Fig.[0 

The techniques are: 

1 . Narrowband on-line filter at redshifted Lya and one broad- 
band filter with the same central wavelength. While neither 
of these observations samples the line or continuum indi- 
vidually, they both sample the line-flux plus continuum; 
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Fig. 1. Schematic diagram of filter positioning for the four 
observational approaches described in Sect. 12.31 the text. 
Narrowband filters (on-line, black) have widths of 2% their 
central wavelengths and broadband filters (continuum, grey) 
have widths of 20% A c . Wavelength indications represent cen- 
tral wavelength of the filters and are marked in the restframe of 
the SED. I.e. observation-frame filters are centred at (1 + z) ■ A c . 
[See the electronic journal article for the colour version of this 
figure.} 

the flux in each filter is the sum of the line-flux (f\, y «) and 
the product of the continuum flux-density {f ont) and filter- 
width (W). If a flat continuum (J3 = 0) is assumed across 
the broadband filter, Fi,y a and f cont can be obtained from 
the solution of two simultaneous equations as 



W h F n - W n F h 
W h -W n 



F h - F n 
W b - W n 



(3) 



where F and W represent the flux and filter width, and sub- 
scripts n and b refer to the broad and narrowband filters, 
respectively. This is the same observational setup as used 



by Cowie & Hu and Fynbo et al. l lIUOll . and divid- 

ing the first part of Equ.[3]by the second, one obtains the 
equivalent width expression used by Malhotra & Rhoads 

2. Narrowband on-line filter centred on redshifted Lya and 
one completely off-line broadband filter centred at (1 + z) ■ 
1500 A to sample the continuum. At 1500A, a 20% (300 A 
full width) filter will not transmit the Lya line itself. The 
continuum flux at Lya can then be estimated by assuming 
the continuum can be approximated by a power-law of the 
form fx oc A 13 , with arbitrary index (3. This technique al- 
lows the investigation of the reliability of assuming the be- 
haviour of the UV continuum slope in the situation where 
only one continuum observation is available. 

3. Narrowband on-line filter centred at redshifted Lya and two 
broadband continuum filters centred at (1 + z) ■ 1500 and 
(1 + z) ■ 2200A. This method assumes the UV continuum 
can be approximated by a power-law of the same form as 
in technique #2. From the two continuum fluxes, the index 
f3 is measured between restframe 1500 A and 2200A. This f3 
is then used to scale the 1500 A flux to the domain sampled 
by the narrowband filter and estimate / cont at 1216A. 

4. Narrowband on-line filter centred at redshifted Lya and 
fo ur broadband filters centred at restframe 1500A, 2200A, 
3300A, and 4400A. These wavelengths are selected to cor- 
respond to the central wavelengths of the filters used in our 
HST/ACS imaging survey of local starbursts (Hayes et al. 
2005 1. In this study we found we could find non-degenerate 
fits in age - Eb-v space if we sampled the UV continuum 
slope and 4000A break. The same technique is adopted 
whereby we fit the synthetic SB99 spectra to the 'real' data. 
A standard least-squares algorithm is used for different ages 
and internal reddenings using a given reddening law. From 
the best-fitting synthetic spectrum, the flux ratio between 
the (1 + z) ■ 1500 A continuum filter and the narrowband 
on-line filter is computed. This scale-factor is then used to 
scale the observed continuum flux at 1500 A to the contin- 
uum flux in the spectral region sampled by the on-line filter. 

3. The simulations 

Our simulations aim to understand the way in which astrophys- 
ical conditions in the host-galaxy and inter-galactic medium 
(IGM) manifest themselves in the determination of F^ya and 
W^ya. To this end, we devised a number of different tests, vary- 
ing individually the parameters that describe the host galaxies, 
in order to examine their impact. In order to standardise the pa- 
rameter setup, a 'standard' restframe template was defined, the 
key parameters of which can be seen in Tab.^ 

By default, single stellar populations are treated, with the 
inclusion of stellar and nebular emission, the relative contri- 
butions of which are left unchanged. Unless otherwise stated, 
the equivalent width of the Lya line added to the spectrum is 
100A for ease of visualisation (I.e. fractional/percentage devi- 
ations are easily interpreted). No internal or Galactic extinction 
is applied to the spectra by default. 

We determined it necessary to run tests to study the impact 
of redshift and H 1 absorption along the line-of-sight, internal 
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Table 1. The parameters of the 'standard' restframe SED 



Parameter 


Value 


w Lya 


100 A 


Metallicity 


0.001 


IMF a 


-2.35 


IMF M low , M up 


1.0M o , 100M g 


Age 


4Myr 


Burst mode 


Instantaneous 


Composition 


Stellar + Nebular 


Secondary population 


No 



Note. IMF of the form: dN/dM cc M- a dM 



dust reddening, and mixing of stellar populations. Burst age is 
not considered in this paper for the following reason: our pre- 
liminary tests showed that the recovery of Fj_y a and Wiya by all 
four techniques were self-consistent to within 2%, (no age de- 
pendence on the recovered observables with age) over the first 
200Myr. In addition, Chariot & Fall (119931 demonstrated Lya 
in emission can only be expected from considerably younger 
star-forming galaxies of ages < 40Myr; over this time we found 
evolution in the observables below the 1 % level. 

Technique #2 allows for the UV continuum slope (J3) to 
be arbitrarily chosen. In a large sample of LBGs re-sampled 
into quartiles according to W]_ ya , Shapley et al. ( 2003 1 find me- 
dian values of j3 = - 1 .09 for the strongest Lya emitters (upper 
quartile) with the continuum flattening to /3 = -0.73 for the 
weakest, lower quartile where the median WLya is -14.92A. In 
contrast, the Starburst99 models show /3 in the range -2.6 to 
-2.0 over the first ~ 50Myr. Excess flux in a narrowband fil- 
ter (n - b < X) is a typical selection function although the the 
amount of excess X is arbitrarily chosen. X = would cor- 
respond to a flat continuum, whereas X < would imply the 
continuum level is rising towards Lya. While X can be com- 
puted by making some simple assumptions, the range in the 
power-law index (J3) from which X would be computed can be 
very large. Without any additional data to constrain the con- 
tinuum slope, any assumption is ad hoc but, for the sake of 
simplicity, in the simulations presented here we assume /3 = 
(i.e. a flat continuum). Moreover, only a modest dust content 
is required to flatten the steep continuum slopes observed in 
the FUV. Additionally, it should be mentioned that the Equ.[3] 
(technique #1) only holds exactly if /? = 0, although the total 
spectral region sampled by this technique is much smaller. 

3.1. Cosmic distance and the IGM 

Any two restframe SEDs generated from the same initial pa- 
rameters are identical. The only way two identically created but 
redshifted SEDs can differ is via the differing effects that result 
from intervening H i in the IGM. The z-distribution and column 
density of H i clouds is pseudo-randomly generated, and as a 
result the redshifted SEDs can differ greatly. Thus a Monte- 
Carlo-type approach is adopted. A population of 10000 rest- 
frame SEDs is generated which, after redshifting, differ only 
bluewards of Lya (except for the possible red wing of DLA 



systems that may extend on to the red side of the line - the 
absorption line-centre is still on the blue side). For this pop- 
ulation, mean and median averages of the fluxes and equiva- 
lent widths are then computed using the techniques described 
in Sect. 12. 31 These computations are performed for populations 
at a range of redshifts between zero and z - 9.0, the highest 
redshift at which narrowband survey techniques have been em- 
ployed in an attempt to find LAEs (Willis & Courbin, 2005). In 
a real imaging survey, there must be additional selection crite- 
ria by which a candidate is determined to be an emitter. In this 
case, we also recompute these mean values after the rejection 
of all objects for which restframe Wi, ya value was determined 
to be less than 20A. This value is hereafter referred to as the 
"refined mean". 

It would also be interesting to simulate more realistic sam- 
ples of galaxies. The Lya luminosity function has been com- 
piled at z = 5.7 and 6.5 by Malhotra & Rhoads H2004L 
However, while spectroscopic observations were used in the 
compilation of these luminosity functions for statistical pur- 
poses, the luminosities of the LAEs themselves have been 
derived photometrically. Hence, they may well contain ob- 
servational biases of the type we attempt to address here. 
Additionally, we would have to assume the characteristics of 
the observation itself such as sky-noise and it was deemed more 
revealing to study populations of identical objects. 

3.2. Internal dust reddening 

In the absence of intervening H i absorption, we study the ef- 
fect of internal dust extinction on the reliability of our recovery 
of Fi,y a and W^ya- Using the common starburst extinction laws 
[ Starburst: Calzetti et al. (2000 1; SMC: Prevot et al. (IT985V 
LMC: Fitzpatrick et al. ( 1985 1 ] we artificially redden the tem- 
plate SEDs and investigate how reliably the various techniques 
estimate Lya and W^ya- Technique #4 relies upon fitting age 
and internal reddening in order to estimate the continuum flux 
at Lya. In these experiments, we also investigate the reliabil- 
ity of this technique when different extinction laws are used for 
the SED fitting from the one used to redden the intrinsic spec- 
trum. This allows us to investigate how well Lya observables 
are recovered if an inappropriate reddening law were assumed. 
For example, we redden the intrinsic SED with the SMC law 
and use the Calzetti law in the SED fitting routines. We then 
investigate whether possible improvements can be made by in- 
cluding the reddening law itself as a free parameter in the fitting 
routine. Tests are carried out in the range 0.0 < Eb-v < 0.5 at 
redshifts of 2, 4, and 6. 

3.3. Underlying stellar populations 

Again without the effects of intervening Hi systems, we in- 
vestigate how accurately and W^ya are recovered when 
the Lya emitting object hosts varying components from an un- 
derlying population. This is mainly relevant to technique #4 
which relies upon sampling of the 4000A break, and, while we 
don't study Fi, ya , or W]_ ya vs. age directly, this study addresses 
the effect an aged stellar population may have. As before we 
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adopt the default template spectrum as a 'primary' SED (I.e. 
the single stellar population starburst spectrum defined in Tab. 
^with an age of 4Myr). We then add an assortment of older, 
post-starburst spectra to the primary, scaling to varying nor- 
malisation factors at a wavelength of 4500A (i.e. scaling by the 
approximate B-band luminosity). We thereby artificially cre- 
ate young starburst galaxy spectra with a population of older 
stars. This normalisation factor is designated «45oo and, by this 
definition: «45oo = is just the default SED defined in Tab.^ 
"4500 = 1 means the default and aged populations contribute 
equally at 4500A; and «45oo = 10 means the old population is 
a factor of 10 more luminous than the starburst SED at 4500A. 
The parameters used in the generation of SED of the under- 
lying population can be seen in Tab. [5]but briefly, we perform 
tests to examine the age of the underlying population, its metal- 
licity, and its normalisation coefficient, «45oo- Then, after red- 
shifting the SEDs, we investigate what values of Fiy a and W^ya 
our techniques recover. When the SED-fitting technique is used 
(technique #4), we always used the unmodified set of starburst 
templates for fitting. Tests are performed at redshifts of 2, 4, 
and 6. Tab. [2] shows the different old stellar populations and 
their contributions as applied to the standard setup, along with 
the resulting W]_ ya calculations. 



4. Results and discussion 

4.1. Cosmic distance and IGM 

As described in Sect. l2.3l and Sect. l3.ll we generate populations 
of galaxies at a given redshift, and determine the mean, me- 
dian, and refined mean values of observed Fiy a and W^ya for 
this population using the four observational techniques. Fig. 
|3 shows a set of histograms representing the distribution of 
these observed quantities for a population of 10000 galaxies 
at z — 5.7, selected to correspond to the redshift of number 
of Lya surveys (eg. Thommes et al. 1998; Rhoads & Malhotra 
I200T1 Westra et al. EBDSl . Restframe W Lya is 100 A. 

A noteworthy feature, present in all these plots is the pop- 
ulation of ~ 750 galaxies with Fi, ya and W^ya around zero 
(and slightly negative). In the line-of-sight to these objects a 
DLA system has been generated near to the target LAE, the 
red damping wing of which has removed the Lyo- line from the 
spectrum. The slight negativity of these values comes from the 
fact that, in order to estimate the line-only flux, the continuum 
has been subtracted and, in all these cases, continuum subtrac- 
tion has resulted in negative fluxes in the line. These objects 
would not be found by imaging observations that target emit- 
ters and hence including them in the computation of average 
quantities would be misleading. This effect therefore leads to 
a redshift-dependent incompleteness, which at z = 6 is about 
10% for these filters. Such effects would need to be treated in 
the computation of the Lya luminosity function. It is for this 
reason that we compute refined mean values for the population 
of objects with detected WLy > 20A only. 

Dashed vertical lines have been added to all these his- 
tograms in three colours. The red line shows the flux or equiva- 
lent width of the Lya line that was added to the restframe spec- 



trum, the green line shows the mean of the recovered quantities, 
and the blue line shows the refined mean. 

The top pair of plots in Fig. |2] represent the quantities as 
determined using technique #1 (a narrow and broad filter cen- 
tred at (1 + z) ■ Aiya). The F^a histogram shows a strong peak 
in the distribution that perfectly corresponds to intrinsic line 
flux. The mean of the distribution falls 15% short of this value. 
The coincidence of the peak in flux distribution and intrinsic 
flux demonstrates the power of this method in determining line 
flux. The top right panel in this figure shows the distribution 
of W]_ ya when computed by this technique (#1). While F^a 
is accurately recovered and evenly distributed around the cor- 
rect value, the Wt, ya distribution cuts on sharply at WL yff =100A, 
peaks, and exhibits a long tail, extending to 220A. While some 
flux may be removed from the narrowband filter by intervening 
H 1 clouds, the broadband filter samples the continuum signif- 
icantly bluewards of Lya. Hence H 1 clouds in the IGM can 
remove a significant amount of flux from this broadband obser- 
vation. Consequently, the continuum flux at Lya can be drasti- 
cally underestimated, resulting in a doubling of observed W\_ ya 
in extreme cases. The black dotted line in this plot shows the 
fraction of objects with recovered Wi, y „ greater than the value 
on the abscissa. It crosses the red dashed line (100A) at 0.8, 
implying that W\_ ya has been overestimated for 80% of all the 
objects (including those with WL ya ~ that represent the in- 
completeness in the observed population). After removing the 
Wi, ya ~ OA objects, this fraction is more like 90%. 

The second pair of plots represents the distribution of the 
same detection properties, as determined by technique #2. Here 
the same on-line filter is used but a single broadband filter cen- 
tred at ( 1 +z)- 1500A is used to estimate the continuum at Lya by 
assuming /3 — between the spectral regions sampled by these 
filters. Since our standard restframe template is that of a 4Myr 
starburst with no dust, the continuum is increasing towards the 
line (J3 ~ -2.6) and the subtraction of an underestimated con- 
tinuum results in a slight shift of the distribution towards higher 
fluxes. In this case where the line is strong, the measured flux in 
the narrowband filter is dominated by the line, and the relative 
effect of underestimating the continuum is small. We have veri- 
fied that the effect is more significant when smaller rest equiva- 
lent widths are input (i.e. continuum-dominated observations). 
In this example, underestimating the continuum flux has a more 
pronounced impact upon the estimate of W]_ ya and the mode in 
the W Lya distribution is shifted from 100A to ~ 170A. While 
technique #2 may provide reliable line-fluxes when the line 
is strong, any equivalent widths derived in this manner can- 
not be considered robust. Even if there is no line {W]_ ya — 0), 
for a young, unreddened burst with f3 < 0, the assumption of 
/3 = will result in a positive observed equivalent width and 
false positive detection of narrowband excess. 

Techniques #3 and #4 reliably reproduce both F]_ ya and 
W]_ ya in this study. Technique #3 assumes the UV continuum 
slope to be a power-law between 2200A and Lya and uses off- 
line filters placed at restframe 1500 and 2200A to extrapolate 
to Lya. Technique #4 makes no assumptions about the con- 
tinuum, and uses a 4-filter SED-fitting technique to estimate 
the continuum level at Lya. With the exception of the small 
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Fig. 2. Histograms showing the distribution of detected Lya flux (left) and equivalent width (right) for a population of 10000 
galaxies at a redshift of z — 5.7. The descending pairs of plots represent distributions as computed by observational techniques 
#1, #2, #3, #4. The vertical dashed red lines indicate the true line flux and equivalent width (I.e. that added to the restframe SED). 
The vertical green line shows the mean values computed for the whole sample of 10000 objects, while the blue line shows the 
mean value of these quantities for all objects with Whya> 20A (refined mean). The dotted black line in the uppermost equivalent 
width plot shows the fraction of galaxies in the sample with observed Wh ya greater than the value on the abscissa axis (I.e. for that 
value of Wl yu ., the fraction of galaxies that have W]_ ya greater than this value). [See the electronic journal article for the colour 
version of this figure.] 



peaks around Fjjy a and Wt, ya -0, these distributions are strongly 
peaked around their intrinsic values. The narrow distribution of 
W^ya demonstrates how successful these techniques are at es- 
timating the continuum fluxes at line-centre. The refined mean 
values of W]_ ya are in error by only ~ 5%. 

In reality, all the accuracy of all these results would be de- 
pendent on the recovered flux in the various bands and the sky- 



noise. While this is not strictly the angle we have chosen for 
this article, we investigate and discuss the effects of sky noise 
below. Fig.|4]demonstrates how the various technique fare near 
the detection limit with the the application of a simple sky- 
noise model. 

Fig.|3shows the evolution of the total mean, refined mean, 
and median values of F Lyo and W^yo as a function of redshift. 
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The blue curves represent technique #1, green curves technique 
#2, red curves technique #3, and black curves technique #4. 
The flux values are the average observed fluxes, normalised by 
the computed intrinsic line flux at the observatory (i.e. the in- 
trinsic line-luminosity scaled down by the luminosity distance). 

As previously discussed, technique #1 reliably recovers 
Fiy a in the line-dominated limit and this is here demonstrated 
to be independent of redshift. However, Wiy a as computed us- 
ing this technique becomes progressively worse with increas- 
ing redshift as can be seen from the right plot - the refined 
mean value for W^yo atz = 6.5 is overestimated by 20%. This is 
a result of the increasing redshift-density of intervening clouds 
with redshift (oc (1 + z) 2 ' 45 ), causing a systematic decrease of 
the estimate of the continuum flux. Technique #2 overestimates 
Fhya by ~ 10% at z = 2. This results from the assumption of 
a flat spectrum not accounting for the blue continuum slope of 
a young stellar population. This technique is unreliable in esti- 
mating Wtya and at all redshifts as the right plot shows, never 
comes close to estimating the correct W^ya- Again as demon- 
strated in Fig. 13 techniques #3 and #4 both reliably determine 
Fiya and W]_ ya at low and high redshift - the refined mean val- 
ues of W]_ ya deviate from the expected values by less than 9, 
and 5%, respectively, at z — 6. 

It is, of course, possible to fit a function to any of 
the redshift-evolution curves shown in Fig. |5J thus obtain- 
ing a functional prescription for the fractional under- or over- 
estimate of F^ya and W^ya as a function of redshift. Such a for- 
mula would be highly desirable as it would be directly applica- 
ble to the results obtained by previous surveys. However such 
a corrective formula would not only be a function of redshift, 
but also of the filter widths, response profiles, and positioning 
in wavelength. Of course, such information is well-known and 
we recommend simulations to be carried out on a survey-by- 
survey basis, and with appropriate consideration of the errors 
and selection criteria. 

Malhotra & Rhoads (2002i used a technique similar to 
technique #1 to observe LAEs at a redshift of 4.5, uncovering 
a population of galaxies with curiously high equivalent width 
(median W]_ ya ~ 400A). Since there is no AGN activity associ- 
ated with these objects (Wang et al. 2004 1 and the maximum 
Wi,ya available from star-formation (with 'normal' IMFs and 
metallicities) is 240 A (Chariot & Fall 1993 1, such a discovery 
attracts speculation. Using our simulations we determine that at 
z = 4.5, the median value of Wt, ya may be overestimated by a 
factor of only ~ 12% - clearly Malhotra & Rhoads' high-WL y » 
result is not an observational effect of the type under consid- 
eration in this article. However, the effect of intervening H i on 
Wijy a as computed by technique #1 is that a large spread in W^ya 
is produced (see Fig. [2 top right), leading to a significant frac- 
tion of LAEs with overestimated Wt, ya - 20% of all the objects 
have W]_y a overestimated by more than 40%. Certain selection 
functions (eg. narrowband excess) will then be more inclined 
to pick up these objects. 

In order to further investigate some more realistic obser- 
vational effects, we implemented a simple noise model. Taking 
the fluxes obtained from the population of objects shown in Fig. 
12(10000 objects at z = 5.7 with W Lyff =100A), we assumed 



that the criterion for a 'detection' is 5cr. We then assigned a 
S /N = 5 to each narrowband flux, and weighted this by the ra- 
tio of detected F\jy a to the intrinsic line flux (Fl^o)- S /N was 
assigned to the broadband observations in an identical manner: 
5 weighted by the ratio of the observed flux to the intrinsic flux. 
Note this weighting modification of S /N = 5 is only applicable 
to the broadband filter centred at Lya - the redder filters are not 
affected by H i clouds in the IGM, and all have S/N = 5. With 
cr for all the observations, we used Box-Muller transforms to 
generate Gaussian deviates around the fluxes in each passband. 
We fed these fluxes back into the expressions used to derive 
Wis/a and replotted the resulting distributions of W]_ ya which 
can be seen in Fig.|4] 

In all cases, the main distribution is now significantly 
broadened; the mode has shifted towards lower values of W^ya, 
and a high W^ya tail has been produced. This results from 
the fact that equivalent width is the quotient of two values 
(Fhya/fcont)- Symmetrical redistribution of the denominator re- 
sults in asymmetric redistribution of the combined quotient; 
compressed on the low side of the mean (lowering the mode) 
and extended at high values. This is most striking using tech- 
nique #1 since signal in the broad filter is lost to H i absorption 

in the IGM 20% of objects have Wi, ya overestimated by a 

factor of 2. Still the number of objects scattered to very high 
Wi, ya is small overall, we reiterate that some selection crite- 
ria will include them: in a narrowband survey they are likely 
to exhibit clear narrowband excess. Redistribution of Lya has 
caused the plots for all techniques to appear similar: modes 
have been reduced slightly and extended, high-WL yff tails are 
present. One noteworthy feature about the progression from 
technique #1 through #4 is the movement of the refined mean 
of the distribution (blue vertical line) towards the intrinsic value 
of 100A, and for the same reason, the steepening of the black 
dotted line. The addition of more filters is necessary to prevent 
the extreme spreading of the distribution and scattering to ex- 
treme values of Wt, ya - For technique #4 the refined mean value 
now accurately recovers the intrinsic values of W^ya, and Wt, ya 
is overestimated by more than 50% for very few objects. 

4.2. Internal dust 

The tests presented here showed indistinguishable results at 
redshifts of 2, 4, and 6, hence the results we present regard- 
ing reddening can be assumed to be independent of redshift. 
Fig- E shows how well W^ya is recovered for some test cases 
when dust is added to the system. The internal reddening for 
a given extinction law (type of dust) is represented on the ab- 
scissa of the plots, while the recovered value of W^yo is pre- 
sented on the ordinate. Where the SED-fitting technique (#4) is 
concerned, the assumed extinction law used by the SED-fitting 
routine may be different from that used to redden the restframe 
SED and can be seen in the central caption of each figure. In 
the tests presented here we only redden the restframe SEDs 
using the Calzetti and SMC laws because all other laws show 
a 2 175 A graphite feature. If we believe we are dealing with 
objects with a very strong UV radiation field, larger graphite- 
based molecules will be destroyed, removing any 2 175 A fea- 



10 



Hayes & Ostlin: The detection properties of high-z Lya galaxies 




2468 2468 
redshift (z) redshift (z) 



Fig. 3. The evolution of F^ya and Wt, ya as function of redshift. The left plot shows the recovered F^ ya , normalised by the predicted 
flux at the observatory, calculated from the intrinsic line luminosity added in the restframe (I.e. the fractional deviation from the 
true value). The right plot shows the evolution of W]_ ya - Colour blue represents technique #1, green: technique #2, red: technique 
#3 and black: technique #4. Dotted lines represent total mean values, dashed lines represent median values, and solid lines show 
the refined means. [See the electronic journal article for the colour version of this figure.] 
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Fig. 4. Histogram showing the detected W\_ ya of 10000 galaxies at z = 5.7. These are the same objects as plotted in Fig.|2]but 
fluxes have been regenerated within assumed observational errors from a simple assumed noise model. Dotted line and vertical 
dashed lines as in Fig. 13 [See the electronic journal article for the colour version of this figure.] 



ture from the reddening vector and producing a law similar to 
that of the SMC (Mas-Hesse & Kunth 1999). We now only con- 
sider the impact upon the determination of Wh ya , having shown 
it to be a much more sensitive performance indicator than F^ ya . 
Intervening H i absorbing systems have been "switched off' for 
this experiment. The range in 0.0 < Eb-v < 0.5 is chosen. 



Technique #1 (blue line) is here demonstrated to very accu- 
rately recover W\_ ya in all cases shown. At Eb-v = 0.5 this 
technique overestimates W\_. ya by just 3% using the Calzetti 
law (central panel) and the technique is largely insensitive to 
dust. This minor effect is a result of the steepening extinc- 
tion curve removing flux from the broad bandpass, leading to 
a mild overestimate of the continuum flux at Lyo-. Conversely, 
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Fig. 5. Recovered Wt, ya as a function of internal reddening using the starburst extinction laws. The abscissa shows the internal 
reddening Eb-v, using a known extinction law. The ordinate shows the recovered value of W^ya using an an extinction law 
assumed by the observer. Left: True dust: SMC-law; assumed dust: LMC. Centre: True dust: Calzetti-law; extinction law free 
parameter in SED fitting routine Right: As for centre but with noise model applied. Colours represent observational strategies as 
in Fig. |3 The filter setup is as described in section l2~3*l (soIid lines). The dashed lines show the same results but with the nearest 
off-line filter (restframe 1500A) shifted to 1400A. [See the electronic journal article for the colour version of this figure.] 



this technique underestimates Wiy a by 13% for SMC-type dust 
(left panel). Across the wavelength range of the 1216A-centred 
broad bandpass (~ 1000A to ~ 1340A), the SMC-law provides 
more extinction than the Calzetti law for a given Eb-v, and 
has a steeper gradient (more rapidly increasing with decreasing 
wavelength). This results in the stronger suppression of the Lya 
line with the SMC law than that of Calzetti, and turns the mild 
overestimate of W^ya into a more significant underestimate. 

Technique #2 (green lines) is highly ineffectual at repro- 
ducing the intrinsic equivalent width when dust is added to the 
system. This is also dependent upon the chosen extinction law: 
SMC-type dust (left panel) completely suppresses the unred- 
dened excess of a Wi, ya = 100A line relative to the 1500A flux 
by E B -v = 0.48. That is, using technique #2, at E B -v > 0.48 
an LAE of W^y a - 100A will be seen as a Lya absorber. This 
effect is nowhere near as extreme when the Calzetti law is ap- 
plied since the SMC law is significantly steeper in the FUV. 

Modeling the UV continuum (J3) as a simple power-law 
(technique #3; red lines) and extrapolating to Lya only per- 
forms marginally better in the presence of dust. A clear down- 
ward trend is visible with Eb-v, although the dependency is 
highly sensitive to the extinction law. At the relatively modest 
value of Eb-v = 0.1, a Lyo- emitting object will have W^ya un- 
derestimated by around 10% using the Calzetti law but by 25% 
for SMC dust. By Eb-v = 0.5 using the Calzetti law, W^ya is 
underestimated 30%, while for SMC dust the line has been al- 
most completely suppressed. Extrapolating /3 does not provide 
a reliable estimate of W\v a when even a modest amount of dust 
is present since the dust extinction curve modifies /3 in such a 
way that it becomes inconsistent with a power-law approxima- 
tion. The dashed lines in Fig. [5] show how recovery of Whya 
is only very slightly improved when the continuum filter sam- 
pling 1500A is moved to 1400A. Of course, it stands to reason 
that moving the off-line filter nearer to the line is going to yield 
a better result and this is something that must be considered 
when designing any observation - even when the filters can 
be ideally placed (which is also subject to the presence of sky 
lines), the gain may not be significant. 



Clearly technique #1 is not greatly susceptible to the ef- 
fects of internal dust reddening. In contrast, techniques #2 and 
#3 are highly susceptible to dust and these observational con- 
figurations do not provide enough information to handle any 
reddening in a reliable manner. Continuum extrapolation tech- 
niques are not sufficient and a technique is required that either 
handles the reddening explicitly or covers a sufficiently nar- 
row spectral domain. Indeed, the technique #4 also becomes 
increasingly unreliable with increasing dust content if the ex- 
tinction law/dust type is not well known - assuming the wrong 
dust type renders the SED fitting method ineffective (Fig|5] left 
panel). However, the 2200A continuum filter is well positioned 
to sample the 2 175 A graphite feature in the reddening vector. 
This enables a third dimension to be added to the SED-fitting 
routine: the discrete parameter of the dust extinction law. By 
fitting age, extinction law, and Eb-v, we are able to recover 
Wiya to within 1 % for all values of Eb-v for each of the extinc- 
tion laws. See the central panel of Fig.[3]for and example using 
Calzetti law. 

Due to concerns about fitting three parameters with noisy 
data, we adopted a similar noise-model to that described in 
Sect. 14.11 Of particular concern was whether we could ac- 
curately recover the extinction law. All fluxes were assigned 
S IN = 5, randomised within the corresponding Gaussian cr, 
and plugged back into the formulae to compute Wt, ya - We now 
generate 1 000 LAEs for each point in Eb-v, and compute the 
mean values recovered. Since there are no Hi IGM clouds in 
this simulation, there is no non-recovered population to that 
needs to be taken into consideration. The results can be seen in 
the right panel of Fig. [5] Clearly, at the 5<x limit, the concerns 
about recovery of the reliability of the fitting routine are not 
serious - the SED fitting routing recovers W^ya to within 3% in 
all cases. 

Obtaining such observations of continuum-faint objects at 
high-z requires a substantial investment of time; many narrow- 
band Lya surveys find a significant population of objects with 
no apparent UV or optical continuum. This is one of the areas 
in which the class of extremely large optical telescopes has the 
opportunity to make a significant impact: with gains in collect- 
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ing area of a factor of 25, restframe UV detections of dusty 
systems may become more of a possibility. This is one of the 
reasons why we push this study to seemingly extreme values 
of Eg-v- If detections can be made, this will be of importance 
given that many Lyo- blobs (LABs) have been detected by nar- 
rowband imaging observations. A significant fraction of these 
LABs are also bright sub-mm sources (Chapman et al. 2005 
Geach et al. 2005 ) implying a very significant dust content and 
therefore reddening of the stellar continuum, although it has not 
been demonstrated that the Lya and sub-mm radiation originate 
from the same regions of the galaxy. This is particularly true if 
the Lya production in LABs is the result of accretion of cold 
gas onto dark matter halos (eg. Haiman et al. 2000 1 in which 
case the Lya may be emitted over significantly extended areas. 
There may also be a certain degree of decoupling between Lya 
photons and the nearby UV continuum due to resonance scat- 
tering. In the cases where there is no continuum detection at 
all in any bands, (eg. Nilsson et al. 2006 1 then equivalent width 
has no meaning. However, if there is a detection of (potentially 
reddened) stellar continuum then total SED will be the super- 
position of this with the gas spectrum that gives rise to Lya, 
whatever the mechanism for its production. In this case, the 
SED fitting approach should be applicable and Lya equivalent 
widths should be recoverable. 

4.3. Underlying stellar populations 

Again intervening H i absorbing systems have been switched 
off for these experiments and we only consider the Lya equiv- 
alent width. Tests here are carried out using the standard tem- 
plate spectrum (defined in Sect. [3] Tab.^ combined with vari- 
ous old stellar populations at redshifts of 2, 4, and 6. However, 
within the tests performed here, no redshift dependence at all 
was detected and hence, none will be discussed. The results 
presented here can be considered to hold at all redshifts. 

Tab. |2 shows the way in which varying the contribution 
from an old stellar population affects the computed values of 
W^ya when the age and metallicity of the underlying popula- 
tion are varied. Fig. ^demonstrates the effect of increasing the 
relative contribution of a 300Myr underlying population. 

For the tests presented for the age and metallicity of the un- 
derlying population, «45oo is set to 1 in every case, meaning the 
underlying population and starburst SEDs contribute equally 
in the B-band. Clearly technique #1 is entirely unaffected by 
the presence of an aged stellar population, being in error by 
no more than 0.5% in every case. Again, the strength of this 
technique results from the fact it is dependent only on a very 
narrow spectral region. Technique #2 is consistently inaccu- 
rate, showing the typical overestimate of W^ya by ~ 70% but is 
self-consistent to within about 10%. Technique #3 is accurate 
to within 10%. The 4-filter SED-fitting technique (#4) imple- 
mented here is included in this study because we have demon- 
strated it to be a reliable technique by which the continuum 
at Lya can be accurately estimated (Hayes et al. 2005 1. The 
strength of this technique relies on the UV continuum slope 
and 4000A break being adequately sampled: then the effects 
of age and reddening aren't completely degenerate and can be 



Table 2. The influence of old stellar populations on detected 
- the effects of age and metallicity. 



Age 


z 


"4500 




W]_y a from technique 




(Myr) 






#1 


#2 


#3 


#4 






0.0 


99.581 


179.774 


95.819 


100.000 


900 


0.008 


1.0 


99.579 


178.865 


98.398 


98.610 


600 


0.008 


1.0 


99.524 


175.017 


96.486 


96.282 


300 


0.008 


1.0 


99.530 


168.273 


92.709 


92.202 


200 


0.008 


1.0 


99.552 


165.828 


91.587 


90.723 


100 


0.008 


1.0 


99.557 


164.678 


90.934 


90.027 


300 


0.001 


1.0 


99.511 


171.039 


94.156 


93.875 


300 


0.004 


1.0 


99.539 


168.296 


92.505 


92.216 


300 


0.020 


1.0 


99.478 


171.927 


94.899 


94.412 


300 


0.040 


1.0 


99.550 


176.421 


97.581 


97.131 



Note: Normalisation factor (714500) defined in Sect. 13.31 
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Fig. 6. The effect of varying «45oo on the recovered values of 
W]_y a as determined by the various techniques. 714500 is defined 
in Sect. 13.31 The error bars on the black line show the standard 
deviation in recovered Wh ya when 1 000 objects are generated 
using the noise model described in the text. [See the electronic 
journal article for the colour version of this figure.] 



disentangled. Naturally, if the 4000A break is so important, it 
is possible that the presence of an additional stellar population 
that contributes strongly to the 4000A break but not in the UV 
could render this technique ineffective. No real trend of Wi, ya 
with metallicity is observed with any of the techniques. A trend 
of W]jy a with the age of the underlying stellar population does 
emerge when computed using techniques #2, #3, and #4: recov- 
ered Whya decreases as the underlying population gets younger. 
The populations are normalised in the restframe at 4500A but 
the SEDs become bluer with decreasing age. Hence, applying 
a younger underlying population has the effect of making the 
composite SED bluer than when an older population is added. 
As a result, the younger the underlying population, the higher 
the flux at 1500 A which leads to a tendency to overestimate the 
continuum at Lya. Therefore the estimates of the continuum 
level at Lya become higher and the computed equivalent width 
is decreased. 
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When considering dependencies upon the normalisation 
coefficient, no trends emerge when using technique #1, (for the 
same reasons this technique is not sensitive to the age of the 
underlying population; Fig. [6ji. Technique #3, which extrap- 
olates the measured value of B to Lya, tends to overestimate 
Whya with increasing 714500 . As W4500 increases, the flatter SED 
of the old population has the effect of flattening in the spec- 
trum in the 1500-2200 A region. However, the 300Myr burst 
turns over towards Lya much more sharply than the very young 
SED, and hence doesn't contributed much at 1216A. Hence the 
continuum flux at Lya is increasingly underestimated as «45oo 
increases. There is a minor tendency for technique #4 to be- 
come less accurate with increasing 724500, underestimating Wh ya 
by 6% at 714500 = 10. The fitting software selects increasingly 
older template spectra as the old population contributes more 
and more to the 4000A break. However, with increasing W4500, 
the composite spectrum changes faster across the 4000A break 
than it does in the FUV. So as older (redder intrinsically) spec- 
tra are selected by the fitting software, in order to maintain a 
good fit in the UV, the best-fit SED will be less reddened than 
it would otherwise. This leads to a tendency to slightly overesti- 
mate the continuum flux at Lya, and a trend of underestimating 
Whya with increasing contribution from an older population. 

In light of these results we made a number of attempts to 
confuse our SED-fitting software. The first was to, as in pre- 
vious sections, see how the SED fitting routine fared when the 
simple noise model was applied. After applying the aged pop- 
ulation at each value of 724500, we randomised all the computed 
fluxes by assigning S/N = 5 in each filter, and computed W\jy a 
from these values. We generated 1 000 objects and computed 
the mean values although these were never found to be deviant 
from the values computed without the noise model at all 714500- 
The standard deviation of the recovered WL yff is represented by 
the error-bars on the black line (technique #4) of Fig. [6] 

Secondly, we randomly selected old populations from any 
of the templates with ages greater than 300Myr, reddened them 
with random Eb-v in the range 0.0 -1.0, and added them to the 
template SED with varying 714500 (range: 0.0 - 3.0). By apply- 
ing two old populations to our template spectra in this fashion, 
we were not able to make technique #4 reproduce W^ yo that 
was in error by more than 10%. 

The final test devised to confuse the software was to vary 
the contribution of the nebular continuum in the template spec- 
tra. The motivation for this test again being that the SED- 
fitting software is sensitive to the Balmer jump. Moreover, our 
understanding of the ionising photon production from mas- 
sive stars is incomplete. There is a lack of data concern- 
ing low-metallicity stars, no single O type stars can be ob- 
served at such FUV wavelengths, and the ionising contribution 
from population m objects is largely speculative. Additionally, 
Lyman-continuum escape has been observed in very few cases: 
Bergvall et al. ( 12006b at low-z in the case of ESO 350-IG38; and 
Steidel et al. (2001 ) at high-z by stacking composite spectra of 
Lyman Break Galaxies (LBGs). It is the reprocessing of these 
photons that determines the intensity of the nebular component 
relative to that of the stars, hence any Lyman continuum pho- 
tons that escape are not reprocessed in the nebular hydrogen 
emission spectrum (lines or continuum). We found that, even 



when scaling the nebular component (see Sect. l2.1> by a factors 
in the range 0-3, technique #4 recovered W]_ ya with errors no 
greater than 5%. These limits extend between two very extreme 
cases: one, where all Lyman continuum photons escape (which 
would result in no nebular emission component and therefore 
no Lya); and two, the possibility that the estimates of ionising 
continuum in stellar atmosphere models are in error by a fac- 
tor of 3 and all the ionising photons are reprocessed as nebular 
emission. 

4.4. Redshift and observational dependencies 

One of the initial motivating factors for the study of Lya, dat- 
ing back to Partridge & Peebles d!967> . was the fact that it 
could provide a tracer of early star formation observable from 
the ground at redshifts greater than around 2. At z = 2, tech- 
nique #1 would be equivalent to a narrowband filter at 3600A, 
plus a U filter. Here the effect of intervening H 1 systems will 
be lower than at higher z and, at the lowest redshifts where Lya 
is observable from the ground, intervening H 1 affects the de- 
termination of W^ya by technique #1 at the 1% level for our 
chosen filter set. That said, any single measurement by this 
technique cannot be deemed reliable without additional obser- 
vations. Followup spectroscopy could obviously confirm the 
presence of intervening systems and would also yield an accu- 
rate measurement of Wh ya , but observational requirements may 
preclude the possibility of obtaining such data for every target. 
Perhaps cheaper, depending on the design of the survey, would 
be to use two off-line filters, say B and V (restframe 1470A and 
1830A, respectively - where CCDs are typically more sensi- 
tive than in U), in order to estimate the continuum level at Lya 
without sampling bluewards of 1216A. However, the Shapley 
et al. J2UU31 LBG sample find median E B -v of 0.099 for the 
strongest Lya emitters with Eb-v increasing with decreasing 
Whya- At Eb-v = 0.1, technique #3 is already underestimating 
Whya by 25% for SMC dust. The dust-sensitivity of technique 
#3 is independent of redshift and will not be discussed further. 
At z — 2, the 4000A break can still be sampled by an observa- 
tion in the /-band allowing technique #4 to be used in order to 
properly handle the effects of dust reddening. 

As redshift increases, the effect of intervening H 1 on tech- 
nique #1 becomes more pronounced, reaching the 10% level 
by z = 4. Systematic uncertainties on numbers derived by this 
technique only should probably be addressed using a Monte 
Carlo-type approach. At this redshift, the 4000A break can 
still be sampled from the ground using the K s band although 
the high night-sky background makes these observations ex- 
pensive. In 40 hours of integration time on an 8 metre tele- 
scope, a typical /f s -band imaging camera can detect an ob- 
ject of K s = 24 in the AB magnitude system at S /N = 5 1 . 
At z — 4, and assuming the local fi-band luminosity func- 
tion (Jones et al. 12006b evolves parallel to the LF at 1500 A 
between redshifts of 0, (Wyder et al. 2005 1 and 4 (Yoshida et 

1 Computations are performed with the Exposure Time Calculator 
for the Infrared Spectrometer And Array Camera (ISAAC) mounted 
on Very Large Telescope (VLT) at ESO Paranal, in ideal observing 
conditions (seeing of 0.7"; airmass of 1.2). 
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al. 2006 1, this corresponds to a detection limit 0.5 mag fainter 
than M* . According to James Webb Space Telescope (JWST) 
Mission Simulator and NIRCam sensitivity estimates, such ob- 
servations to these depths can be obtained at 2pm in just 1 sec- 
ond. 

At redshifts greater than 5, L-band observations would be 
required in order to sample the 4000A break. While deep ob- 
servations may be obtainable from the ground at shorter wave- 
lengths, the sky -background at wavelengths longer than K s pre- 
clude such observations and current groundbased facilities do 
not approach the required sensitivities. Consequently these ob- 
servations would need to be carried out by mid-infrared tele- 
scopes in more unfriendly and expensive environments: for ex- 
ample the Spitzer Space Telescope or JWST, or a proposed 
MIR telescope at Dome-C, Antarctica. This is somewhat con- 
trary to the motivation for using Lya as a probe for primeval 
star-formation. However, we reiterate that it is only this obser- 
vation that needs to be done from space - the other bands may, 
depending on redshift, more economically be performed from 
the ground. Moreover, in the age of deep multi-band surveys, 
the amount of pre-existing data there is for a field is always 
a consideration. If observing a deep field at a certain wave- 
length, it may well be advantageous to select fields for which 
much data is already in existence - observations at new wave- 
lengths are frequently added to existing datasets (the Chandra 
Deep Field, for example). Hence the majority of the observa- 
tions required for technique #4 (targeting certain redshifts, at 
least) may already be in place, requiring only the additional in- 
frared bands and it would be wise for future Lya surveys to cap- 
italise on the currently existing data. It is also likely that early 
JWST programs will include deep observations of regions with 
existing deep optical data. 

In addition to the observational, there are also further theo- 
retical uncertainties; mainly concerning the validity of the stel- 
lar atmosphere models in the UV. These models have currently 
not been well tested at these wavelengths. See Sect. 14.31 for a 
discussion on this. 

5. Conclusions and summary 

Simulations of observations of high-redshift Lya emitting 
galaxies have been performed. Specifically, the simulations 
have examined how efficiently the intrinsic values of Lya flux 
and equivalent width are recovered by narrowband imaging ob- 
servations. F^ya and Wj_, ya are determined by the four methods 
described in Sect.[3]and the effects of intervening H i clouds, in- 
ternal dust-reddening, and underlying stellar populations, have 
been investigated. In summary: 

- Observing Lya with one off-line continuum filter has been 
shown to be highly ineffectual. The steep UV continuum slope 
(j6) may cause overestimates of Wi, ya by a factor of almost 2, 
while the neglection of dust reddening may cause the reverse 
effect: narrowband imaging observations may interpret bright 
Lya emitters {Wt, ya > 100A) as Lya absorbing systems. 

- Using two filters to estimate the UV continuum slope and ex- 
trapolating to Lya using a power-law improves the situation to 
varying degrees. In the dust-free cases this technique becomes 
a highly competitive method of determining F^ ya and Wiya- 



However, even modest amounts of dust render this technique 
highly ineffectual since dust modifies the continuum in such a 
way that the power-law parameterisation becomes unreliable. 
Eb-v = 0.1 results in 25% underestimate of W^ya- 

- Economic observing strategies that utilise a narrowband and 
single broadband filter with the same central wavelength are far 
less susceptible to dust. However, the continuum flux estima- 
tion in such observations becomes suspect due to the blue half 
of the broadband filter being on the blue side of Lya, and hence 
being susceptible to absorption of the continuum by interven- 
ing Hi clouds. Such a technique can lead to overestimates of 
Wi,ya by factors of > 2. With the addition of a simple sky-noise 
model, this effect becomes more pronounced with decreasing 
S IN. We suggest that in such observations, simulations are per- 
formed to estimate the biases caused by such effects. 

- SED-fitting techniques that observe only redwards of Lya 
are not susceptible to intervening H i absorption or dust red- 
dening, when a single extinction law is considered at all wave- 
lengths. Such techniques need to sample the UV continuum 
slope, 2 175 A dust feature, and 4000 A break. This way reliable 
estimates of the continuum flux at Lya can be made and W]_ ya 
estimates are much better constrained in the presence of noise. 

- None of the techniques we considered here were hugely sen- 
sitive to the presence of underlying old stellar populations or 
randomly mixed populations. Additional populations have no 
detrimental effect on the efficiency of the technique utilising 
two filters centred on Lya. The contribution of such popula- 
tions to the 4000A break cause only minor (< 10%) errors in 
the SED-fitting routine. Similarly increasing the relative contri- 
bution of nebular emission has no significant impact upon the 
estimates of Fi^a and W]_ ya . 

- Independent of technique, we find a redshift-dependent in- 
completeness that results from H i systems along the line-of- 
sight. Such an effect would have to be allowed for in the com- 
pilation of a Lya luminostiy function. The red damping wing 
of DLAs close to the LAE can remove sufficient flux to com- 
pletely suppress the Lya line. At z — 6, one could expect to 
miss around 10% of the sources using our perfect filter set. 

We have shown that to recover FL ya and WL ya from high- 
z LAE systems, a SED fitting method is preferable if obser- 
vations can reach the required depths. Other techniques suffer 
badly from either Lya absorption along the line-of-sight or red- 
dening by internal dust. If such off-band observations are not 
available, the use of a single narrow- and broadband filter with 
the same central wavelength seems to be preferable, although 
such observations should be complemented with some tuned 
statistical study of the effects of intervening H i systems. 
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